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ABSTRACT 



Known antenna array manifold information is exploited to 
provide for fast and accurate channel estimation and 
demodulation on both the forward and reverse links, thereby 
increasing capacity in PCS and cellular CDMA networks 
that use adaptive antenna arrays. On the reverse link, an 
"extended" array manifold is used to assist the demodulator 
in maintaining a dynamic estimate of the spatial signature to 
use for beamforming and coherent demodulation. On the 
forward link, channel estimation is performed in the handset 
to provide a robust solution for beamforming. 

48 Claims, 6 Drawing Sheets 
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MANIFOLD ASSISTED CHANNEL 
ESTIMATION AND DEMODULATION FOR 
CDMA SYSTEMS IN FAST FADING 
ENVIRONMENTS 

The present invention is related to, and hereby claims 
priority to, U.S. provisional patent applications Ser. Nos. 
60/097,340 entitled "Practical Space-Time Radio Method 
For CDMA Communication Capacity Enhancement/' filed 
Aug. 20, 1998; 60/097,544 entitled "Integrated Spatial 
Demodulator For CDMA Systems In Fast Fading 
Enviornments," filed Aug. 21, 1998; and 60/111,437 entitled 
"Manifold Assisted Channel Estimation For CDMA Sys- 
tems In Fast Fading Environment," filed Dec. 8, 1998. 

BACKGROUND 

1. Field of the Invention 

The present invention relates to wireless communication 
systems, and in particular, to using adaptive antenna arrays 
in PCS and cellular CDMA networks for capacity enhance- 
ment. 

2. Description of Related Art 

A standard technique used by commercial wireless phone 
systems for increasing capacity is to divide the service 
region into spatial cells. Instead of using just one base 
station to serve all users in the region, a collection of base 
stations is used to independently service separate spatial 
cells. In such a cellular system, multiple users can reuse the 
same frequency channel without interfering with each other, 
provided the users access the system from different spatial 
cells. The cellular concept, therefore, is a simple type of 
spatial division multiple access (SDMA). Note that through- 
out this description, various acronyms will be used, which 
are listed and defined in the Table below. 



Acronym 


Definition 


AOA 


Angle of arrival 


ARIB 


Association of Radio Industries 


BS 


Base station 


CDMA 


Code division multiple access 


CM 


Constant modulus 


DSP 


Digital signal processing 


EMP 


Extended manifold processor 


FDD 


Frequency division duplexing 


FHT 


Fast Hadamard transformers 


FL 


Forward link; 


MAC 


Mu ltip ly-a nd- accumulate 


MAD 


Manifold assisted demodulator 


MMSE 


Minimum mean -squared error 


MS 


Mobile station 


OTD 


Orthogonal transmit diversity 


PN 


Pseudo-noise 


PSC 


Primary spatial correlator 


RL 


Reverse link 


SDMA 


Spatial division multiple access 


SIR 


Signal-to-interference ratio 


SNR 


Signal-to-noise ratio 


SVD 


Singular value decomposition 


TOA 


Time of arrival 



In the case of digital communication, additional tech- 
niques can be used to increase capacity. One well-known 
method is to use spatial signal processing in code division 
multiple access (CDMA) systems. CDMA is normally a 
spread-spectrum technique that does not limit individual 
signals to narrow frequency channels but spreads the signals 
throughout the frequency spectrum of the entire band. Sig- 
nals sharing the band are distinguished by assigning differ- 
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ent orthogonal digital code sequences or spreading signals to 
each signal. Practical techniques for capacity enhancement 
in CDMA systems using an antenna array are described in 
commonly -owned U.S. patent application Ser.No. 08/929, 

s 638 entitled "PRACTICAL SPACE-TIME RADIO 
METHOD FOR CDMA COMMUNICATION CAPACITY 
ENHANCEMENT" and U.S. Prov. patent application Ser. 
No. 60/071,473, entitled "FADING MITIGATION USING 
WIDE APERTURE ANTENNA ARRAY IN SPREAD 

10 SPECTRUM COMMUNICATION SYSTEMS" and 
60/077,979, entitled "CAPACITY ENHANCEMENT FOR 
W-CDMA SYSTEMS", all of which are incorporated by 
reference in their entirety. 

Spatial signal processing can be used in both the forward 

15 link (base station to mobile station) and reverse link (mobile 
station to base station) of a CDMA communication system 
to provide significant signal-to-noise ratio and capacity 
improvements. In the reverse link, spatial signal processing 
includes estimating a spatial signature (defined herein as the 

20 vector of antenna output signals, including multipath 
components, at a given time due to a transmitted signal at a 
certain location, such as described in "Experimental Studies 
of Spatial Signature Variation at 900 MHz for Smart 
Antenna Systems" by S. S. Jeng, G. Xu, H. P. Lin, and W. 

25 J- Vogel, IEEE Trans, on Antennas and Propagation, vol. 46, 
no. 7, July 1998, pp. 953-962, which is incorporated by 
reference in its entirety) of an IS-95 based CDMA signal to 
determine multipath angle of arrival (AOA) values and 
coefficients. An IS-95 system is described in TTA/EIA/IS- 

30 95- A, "Mobile Station-Base Station Compatibility Standard 
for Dual-Mode Wideband Spread Spectrum Cellular 
Systems", May 1995, which is incorporated by reference in 
its entirety. The reverse link coefficients are then used to 
combine a plurality of antenna output signals (after down- 

35 conversion to base band), i.e. beamforming. Thus, the ability 
to accurately estimate the spatial signature is an important 
objective in CDMA systems. However, with a frequency 
division duplexing (FDD) system, the forward link and the 
reverse link occupy different carrier frequencies or bands, 

40 but overlap in time. This difference between the forward and 
reverse link frequencies reduces the correlation between 
fading of the two links so that the two links have signifi- 
cantly different spatial signatures. Therefore, forward link 
beamforming using the reverse link spatial signature esti- 

45 mate is not possible. However, average AOA is generally 
preserved in FDD systems between the forward and reverse 
links for mobile stations far away from the base station. 

Factors which limit the ability to provide for accurate 
spatial signature estimation include the fading rate (Doppler 

50 spread), angle spread, and delay spread profiles of the 
incoming signals. In particular, fast fading, which is created 
by the combination of multipath components of a signal 
being reflected from various elements ("scatterers") in the 
neighborhood ("scattering zone") of a moving transmitter 

55 with random phases, is a major concern in accurate spatial 
signature estimation. The wireless communication channel 
is assumed to have multiple scattering zones characterizing 
the signal propagation between the base station (BS) and the 
mobile station (MS). "Mobile Cellular Telecommunica- 

60 tions" by W. C. Lee, McGraw-Hill, New York, 1995, which 
is incorporated by reference in its entirety, describes scat- 
tering zones around the mobile station. The main scattering 
zone is in the neighborhood of the MS. Large objects, such 
as buildings and waterways, create other scattering zones. 

65 Metal objects in the vicinity of the BS can cause the signal 
to be reflected and can influence the transmission path. 
However, in most cases, the BS antenna array is located 
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above the nearby scatterers, which are assumed to be less pilot signal (e.g., for training) on the reverse link requires the 
significant. As the fading rate or Doppler spread increases, pilot signal to be low power in order to minimize capacity 
the time available to collect coherent data (integration time) loss in the reverse link. However, a lower power pilot in 
decreases. This problem becomes more severe as cellular coherent demodulation requires a longer integration time to 
systems move from the 800 MHz range to the 1900 MHz 5 assure sufficient reference signal quality, 
range or higher, which can increase the fading or Doppler Also, an unknown or varying signal time of arrival (TOA) 
spread by a factor of two or more. For example, a vehicle requires continuous "time searching" and determination of 
moving at 60 mph can induce a Doppler spread in excess of beamformer coefficients at each time hypothesis. In CDMA 
180 Hz in a 1900 MHz system. In general, spatial signature type systems, the signal timing must be recovered before any 
estimation must take place in a duration that is an order of 10 demodulation can take place. Hence, a search process is 
magnitude shorter than the period of the fading rate. For conducted by a series of hypotheses through which the 
example, if the fading rate is 100 Hz, then spatial signature system is varying the time of the reference correlating 
estimation must take place within 1 msec, sequence and then cross-correlating with the incoming sig- 
Various methods for spatial signature estimation or beam- f 1 (f ^ ™ e tin f re ^ d b * each hypothesis must 
former generation have been proposed. These methods can « be sho rt 10 allow a 1™* search (assuming that recovery of 
be characterized by the level of knowledge of the structure tirain | caimot be d ° ne before determmation of beam- 
of the signal impinging on the antenna array and whether or former coefficient since there might not be sufficient signal- 
not a training sequence is present. Characteristics of the ^ noise rat |° * f int)> ™ US ' any f daptlV , C alg °" thm f * 
temporal, spatial, spectral, or modulation structure of the determine the beamformer coefficients must provide fast 
impinging signal may be known and can be exploited in the M convergence in order to achieve a reasonable overall search 
spatial signature estimation, such as described in "Algebraic iimt ' . , . . 
Methods for Deterministic Blind Beamforming" by A. J. Furthermore, MMSE algorithms operating in "decision- 
Van der Veen, Proceedings of IEEE, vol. 86, no, 10, October, directed" mode, as discussed in "Adaptive Filter Theory by 
1998, pp. 1987-2008, which is incorporated by reference in s - Haykin, referenced above, tend to fail at low signal-to- 
its entirety. Many adaptive algorithms based on minimum 25 interference ratio (SIR) conditions at the system input. In the 
mean-squared error (MMSE) or constant modulus (CM) and presence of fast fading, instances of low SIR are common, 
exploiting temporal or modulation signal structure to per- which ™ cause extreme interference or even a "dropped" 
form estimation of the spatial signature are well known, such cal1 - 

as described in "Algebraic Methods for Deterministic Blind Therefore, a system and method are desired that provide 

Beamforming" by A. J. Van der Veen, referenced above, and 30 channel estimation in both forward and reverse links and in 

in "Space-Time Processing for Wireless Communications" fast fading and low SIR environments. 

by A. J. Paulraj and C. B. Papadias, IEEE Signal Processing SUMMARY 

Magazine, vol. 14, no. 6, November 1997, pp. 49-83 and , , . , . . . . 

«. * J iri tu » u c u „v« d™hv^ tun In accordance with the invention, known antenna array 

"Adaptive Filter Theory by S. Haykin, Prentice -Hall, . . . , ' . , - P J . 

c i , n- xt i 100/ u iu ~t\„uLu ,^ • ^™ 35 manifold mformation is exploited to provide for fast and 

Enelewood Clififs, N.J., 1986, both of which are incorpo- A , , . t . j j *, 1 * u *u *u 

, , ' . ' ' . ' 4 r accurate channel estimation and demodulation in both the 

rated by reference m their entirety. £ , , ,. , tL , • • 

7 J forward and reverse links, thereby increasing capacity m 

A disadvantage of these adaptive algorithms is that they pcs and CDMA networks that use adaptive antenna 

generally do not exploit knowledge of the array spatial arrays Qn thc reverse link ^ an « cxte nded" array manifold is 

information or array manifold and generally require sub- ^ Ufjed tQ assist tfae demodljlator in maintaining a dynamic 

stantial time to converge. estimate of the spatial signature to use for beamforming and 

The array manifold is a collection of array response coherent demodulation. On the forward link, channel esti- 

vectors (where each array response vector is the set of mation is performed in the handset to provide a robust 

antenna output signals at a given point in time due to a solution for beamforming. 

far-field signal with no multipath) created by stepping the 45 a de-spread received signal vector y ( - (i-1 to N) from M 

angle of a narrow-band point source (in two or three- antenna element output terminals corresponding to 1 of N 

dimensional space, under far-field and no multipath possible received symbols is spatially correlated with the 

conditions). The array manifold is a trace in M-dimensional array man ifold matrix C to produce a correlation vector p. 

space, where M is the number of antenna elements in the j^q array ma nifold matrix C is an MxK matrix, where M is 

array. Array manifolds are well known and are described in 5Q the number of antenna elements, and K is the number of 

references such as in "Algebraic Methods for Deterministic angles used to create the array manifold (e.g., K~256), i.e., 

Blind Beamforming" by A. J. Van der Veen, referenced each row of matrix C represents one antenna element of the 

above, array manifold, and each column of matrix C represents one 

A signal arriving at the antenna array in a non- multipath angle in the array manifold. This matrix generally spans the 

situation produces a received vector that is on the array 55 whole M-dimensional space produced by the M-element 

manifold. However, when multipath exists, the received antenna array. 

vector is a linear combination of all arriving multipath wave The resulting correlation vector p represents a correlation 

fronts, which is no longer on the array manifold. The between the de-spread, received vector corresponding to one 

Euclidean distance between the array manifold and the 0 f the N possible symbols with each of the K angles of the 

received vector is a function of multipath level, multipath 60 array. The spatial correlation is performed for each of the N 

angle spread, and interference power and increases with an possible received symbols to produce a KxN matrix P, where 

increase in each of these variables. Interference includes the the i th column of the p matrix contains values with magni- 

sum of thermal noise and other incoming transmissions. tudes corresponding to the level of correlation between y ( 

Another disadvantage of some MMSE adaptive algo- and the array manifold for all K angles. Finding the maxi- 

rithms based on nullifying specific interfere rs is that perfor- 65 mum element of the i'* column of P is equivalent to selecting 

mance is reduced when the number of interference sources the column of the array manifold that is closest in Euclidean 

is large (typical for CDMA). Furthermore, using a dedicated distance to the vector y,-. 
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The P matrix is then used to produce a magnitude and stability, and simplify the computational process. The array 

angle of arrival (AOA) estimate each symbol time, which manifold allows accurate AOA and multipath angle spread 

are accumulated to generate an AOA histogram. After low- estimation using "smoothed" data, as described in U.S. Pat. 

pass filtering (or "smoothing"), this histogram will exhibit App. "PRACTICAL SPACE-TIME RADIO METHOD 

"peaks" in the direction of the main scatterers and a distri- s poR CDMA COMMUNICATION CAPACITY 

bution that follows the angular spread of the transmission ENHANCEMENT*' and U.S. Prov. Pat. Apps. "FADING 

source. The multipath scattering area size (or angle spread) MITIGATION USING WIDE APERTURE ANTENNA 

can be estimated from the variance of the distribution around ARRAY IN SPREAD SPECTRUM COMMUNICATION 

a peak. Thus, the histogram can be used to determine the SYSTEMS" and "CAPACITY ENHANCEMENT FOR 

AOA and angular distribution of the most probable signal J(J W . CDMA SYSTEMS", referenced above. The AOA and 

P ■ .... angle spread information can provide "bounds" for spatial 

Next, if the incoming signal angular distribution is signature estimation, thereby eliminating large beamformer 

known, a subspace V spanned by the columns of the array emTS eyen when ^ SIR becomes low . 

manifold matrix C associated with the given angular distri- „ „ , . . „ „. . , . 

bution (mean and variance of the signal AOA) can be „ 0veraU > *f P«sent invention oilers 1) increased imple- 

defined Ttiese columns (vectors) span a subspace that can 15 mentation effic,ency smce the system can ^bc implemented 

be approximated by finding a set ofL orthogonal vectors that cost-effectively within an ASIC not only for base stations, 

. . t f( , . , r . ° u . A . , but also for handsets, 2) high speed operation suitable tor 

contains most of the signal energy. Finding this reduced rank , it _ , • j r v J j i\ . i_i 

u & u j u * i 1 a ™ •** / both stationary and mobile applications, and 3) stable 

subspace V can be done by singular value decomposition (or _ / t lL .i t , . 

a similar procedure) of the original subspace and selecting Performance, i.e., not ^susceptible to convergence and sta- 

the singular vectors corresponding to the larger singular 20 b *ty issues associated with many adaptive spatial filtering 

values. algorithms. 

Since the de-spread, received antenna vector y,. contains ^ P resent invention will be more fully understood upon 

contributions from all directions, it contains energy in the consideration of the detailed description below, taken 

full space spanned by the columns of the antenna array ^ to S eth « r with the accompanying drawings, 

manifold C. Most interference can be filtered out of the BRIEF DESCRIPTION OF THE DRAWINGS 
received vector by projecting this vector y i into the confined 

subspace V. Thus, the sampled received antenna vector FIG. 1 is a partial block diagram of a beamforming system 

(after de-spreading) is projected into this approximated according to one embodiment of the present invention; 

subspace, V, by a simple dot product of the received vector 30 FIG. 2 is a block diagram of a single finger of a manifold 

with each of the columns of matrix V. The projection assisted demodulator (MAD) for the reverse link of an IS- 95 

coefficient vectors corresponding to preliminary data sym- system; 

bol decisions are averaged over a predetermined number of FIG 3 ^ a block diagram of a single finger of a MAD for 

symbols. The number of vectors to average is determined by the reyerse link of 3(J proposed standards; 

the Doppler spread and is selected to correspond to a time « ^ r ^, A . J( c ^.lj j- ui 1 

•a/- u-twu i • * • i i * ut FIG. 4 is a diagram of one of the de-spreading blocks of 

period during which the spatial signature is relatively stable, p[G 3 { ^ ^ w _c DM A system* 

Generally, this period is an order of magnitude shorter than ' y ' 

the period of the Doppler spread. The averaged projection FIG - 5 1S a diagram of one of the de-spreadmg blocks of 

coefficient vector, Z avgl is then multiplied by the approxi- FIG - 3 for a CDMA2000 system; 

mated subspace vectors from matrix V to form the beam- 4Q FIG. 6 is a block diagram of a single finger of a MAD for 

former coefficient vector w. If more than one distinguishable the forward link of a CDMA2000 system; and 

AOA exists, several beamformers are used to track the FIG. 7 is a diagram of one embodiment of the 

multiple scattering zones. de-spreading block of FIG. 6. 

The processed vector w contains less interference and Use of the same reference symbols in different figures 

noise than the original vector y because the subspace pro- 45 indicates similar or identical items, 
jection rejects all the components that are orthogonal to the 

selected subspace. Projecting the received vectors into a DETAILED DESCRIPTION OF THE 

subspace of dimension L reduces the noise power by a factor PREFERRED EMBODIMENTS 

of M/L, where M is the number of antenna elements. If the piG. 1 shows a partial block diagram of a beamforming 

subspace is selected correctly, then the signal power is only 50 system in accordance with the present invention. Through- 

slightly reduced. Thus, signal-to-noise ratio (SNR) for the out t b e description, a beamforming system or demodulator 

beamformer coefficient vector is improved by approxi- implementation that uses the extended array manifold 

mately M/L. This reduces the amount of data that needs to concept, discussed in detail below, is referred to as a 

be integrated for the purpose of estimating the beamformer manifold assisted demodulator (MAD). In FIG. 1, a portion 

coefficients. As a result, performance is improved (e.g., 55 0 f a MAD 100 is shown, which includes a primary spatial 

lower symbol error rate for a given SNR) in fast fading correlator (PSC) 110, a digital signal processing (DSP) unit 

environments when compared with direct averaging of 120, an extended manifold processor (EMP) 130, and a 

de-spread received vectors and compared to conventional beamformer unit 140. De-spread received signal vectors y ( 

two antenna systems with diversity combining. (i=i to N, where N is a selected integer representing the 

The extended manifold concept also can be used to 60 number of possible received symbols) are transmitted to 

provide an initial beamformer coefficient vector for other PSC 110, EMP 130, and beamformer 140. Each signal 

methods such as MMSE adaptive algorithms. Since the vector y,- is a set of M signal components, where M is the 

initial coefficient vector is closer to the desired solution and number of antenna elements in the antenna. Signal vectors 

has better SNR, an adaptive algorithm will converge more y. will be discussed in detail below for IS -95 and "3rd 

quickly. 65 Generation" or "3G" systems (e.g., CDMA2000, as 

Even partial knowledge of the array manifold can signifi- described in "The cdma2000 ITU-R RTT Candidate 

cantly reduce required integration times, provide improved Submission", HA Subcommittee TR-45.5, June 1998 and 
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W-CDMA, as described in "Japan's Proposal for Candidate 
Radio Transmission Technology on IMT-2000: W-CDMA", 
Association of Radio Industries and Businesses (ARIB) 
IMT-2000 Study Committee, June 1998, both of which are 
incorporated by reference in their entirety). Thus, each of 
PSC 110, EMP 130, and beamformer 140 receive N sets of 
M signals y, to y^. The N sets correspond to the results of 
correlation of the de -spread M-element received antenna 
vector against N possible received symbols. FIG. 1 shows 
one "finger" of a Rake- type receiver; the single DSP unit 
120 is shared across all fingers. Rake-type receivers are 
known and described in references such as "CDMA, Prin- 
ciples of Spread Spectrum Communication" by A. J. Viterbi, 
Addison-Wesley, Reading, Mass., 1995, which is incorpo- 
rated by reference in its entirety. Note that PSC 110, EMP 
130, and/or beamformer 140 could also be shared across 
fingers, although the description is for a dedicated PSC, 
EMP, and beamformer per finger. 

Each PSC 110 performs a spatial correlation between the 
MxN matrix Y of de-spread received vectors y. and an array 
manifold matrix C, as given by equation (1). 



p=c H % 



(i) 



where the superscript H denotes a Hermitian-transpose. The 
array manifold matrix C is an MxK matrix, where K is the 
number of angles used to create the array manifold (e.g., 
K=256), i.e., each row of matrix C represents one antenna 
element of the array manifold, and each column of matrix C 
represents one angle in the array manifold. This matrix 
generally spans the whole M-dimensional space produced 
by the M-element antenna array. The \ th column of the KxN 
matrix P contains values with magnitudes corresponding to 
the level of correlation between y,- and the array manifold for 
all K angles. Finding the maximum element of the I th 
column of p is equivalent to selecting the column of the 
array manifold that is closest in Euclidean distance to the 
vector y { . 

Since the correlation process is linear, both relative ampli- 
tude and phase of the incoming signal are preserved. Hence, 
this process can be utilized for both non-coherent demodu- 
lation (e.g., IS-95 M-ary) and coherent phase demodulation 
(e.g., PSK) schemes, both of which are described in 
"CDMA, Principles of Spread Spectrum Communication" 
by A, J. Viterbi, referenced above. 

The matrix P of the PSC 110 is further processed in the 
PSC to produce a magnitude and AO A estimate each symbol 
time as described below for IS-95 and 3G systems. These 
values are sent to DSP unit 120 for generating an AOA 
histogram. Since noise and multipath provide a dynamic 
received vector, a continuous accumulation of AOA samples 
allows an AOA histogram to build up. After lowpass filtering 
(or "smoothing"), this histogram will exhibit "peaks" in the 
direction of the main scatterers and a distribution that 
follows the angular spread of the transmission source. The 
multipath scattering area size (or angle spread) can be 
estimated from the variance of the distribution around a 
peak. Thus, the histogram can be used to determine the AOA 
and angular distribution of the most probable signal paths. 
AOA histogram processing is described in U.S. patent 
application Ser. No. 08/929,638, entitled "PRACTICAL 
SPACE-TIME RADIO METHOD FOR CDMA COMMU- 
NICATION CAPACITY ENHANCEMENT", referenced 
above. 

After creating and processing the AOA histogram, DSP 
unit 120 generates a matrix V, which is an approximation of 
the subspace matrix V. Matrix V is obtained from a singular 



15 
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value decomposition (SVD) of a matrix C 6 , as defined in 
equation (2) below. 

Cq-USV" (2) 

Matrix C 6 has rows that are the columns of the array 
manifold matrix C for the angles Q±o s , where o*/ is the 
angle spread estimate variance (as determined from the AOA 
histogram). Matrix U consists of the left singular vectors of 
matrix C e , matrix S is a diagonal matrix whose entries are 
the non-negative singular values of C e , and matrix V con- 
sists of the right singular vectors of matrix C e . The V 
subspace is approximated by selecting L columns of matrix 
V which correspond to the dominant singular values and 
forming a new matrix V, Dominant singular values are 
defined as those having the largest values. For example, L 
can be selected to include singular values with 90% of the 
total sum of all the singular values. The value L can range 
from 1 to all of the columns of V. Note that in the situation 
where L includes all the columns of matrix V, matrix V is 
equal to matrix V. 

Matrix V is calculated for each distinguishable AOA and 
assigned to the appropriate finger. Matrix V for a particular 
finger is transferred from DSP unit 120 to EMP 130. EMP 
130 then performs a projection of the sampled received 
antenna vector y, (after de-spreading) into the approximated 
subspace, V, by a simple dot product of the received vector 
with each of the S? columns. The projection coefficient 
vector z,. is then given by: 

where y t - (i=l to N) is the Mxl de -spread received vector 
30 corresponding to 1 of N possible received symbols. During 
each symbol time, N possible y ( vectors are projected into 
the subspace V, and the resulting one of vectors z s with the 
maximum energy (defined as HzJ] 2 ) is selected by EMP 130 
for further processing. EMP 130 then averages the selected 
3S Lxl projection coefficient vector z max over a predetermined 
number of multiple symbols (to create z avs ) to further reduce 
the noise. Vector z avg represents the average grojection 
coefficients that when multiplied by the subspace V produce 
an estimate of the spatial signature. Averaging can be 
40 performed on both pilot and data signals. In the case of data, 
the selected vector z max corresponds to a preliminary data 
symbol decision (other techniques for generating a prelimi- 
nary data symbol decision to specify which z i to average 
could also be used). The number of vectors to average is 
45 determined by the Doppler spread and is selected to corre- 
spond to a time period during which the spatial signature is 
relatively stable. Generally, this period is an order of mag- 
nitude shorter than the period of the Doppler spread. Finally, 
EMP 130 constructs a beamforming coefficient vector w by 
multiplying the averaged projection coefficient vector, z , 
by the approximated subspace V, as shown by equation (4). 

H-t/z.^ (4) 

EMP 130 transmits the beamformer coefficients w to 
55 beamformer 140 for processing with the received signal 
response vectors y,. Beamformer 140 performs an inner 
product operation w^y ( .. The N resulting values give an 
indication of the likelihood of each of the N possible 
symbols and are sent to standard decoder blocks as 
60 described in "CDMA, Principles of Spread Spectrum Com- 
munication" by A. J. Viterbi, referenced above. 

FIG. 2 shows a single finger (demodulation channel) of 
the MAD discussed above for the reverse link of an IS-95 
system. In particular, FIG. 2 provides additional detail as to 
65 the processing of the received signal vectors into the MAD. 
The full MAD implementation includes a plurality of MAD 
fingers (at least two, to allow for minimum time diversity). 



50 
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The MAD finger described can perform both time search 
and demodulation. Note that a local buffer memory 200, 
which stores the array manifold data, is shared across all 
fingers, as is DSP unit 120 described above. PSC 110 also 
can be shared across all fingers in an alternative implemen- 5 
tation that would trade fast search time for less complexity. 
Referring to FIG. 2, each of M antenna data streams 
(coming from receiver analog -to -digital converters) is 
de-spread and accumulated, such as described in "CDMA, 
Principles of Spread Spectrum Communication" by A. J. 10 
Viterbi, referenced above. I and Q components of the 
received signal are first fed from an antenna array with M 
elements. The M antenna element output signals are down- 
converted to a baseband frequency and digitized. The M 
antenna streams are then de-spread along M parallel corre- 15 
lation channels as each antenna stream is first input to a 
complex multiplier 210 for multiplication by the appropriate 
long code and then input to a second complex multiplier 220 
for multiplication by an appropriate short code, where the 
long and short codes are supplied from a code generator 230. 20 
A sequencer 235 aligas the timing of the de-spreading 
operation on Walsh symbol and antenna boundaries (antenna 
data may be received serially). The de-spreading process is 
also described in U.S. patent application Ser. No. 08/929, 
638, referenced above. After de-spreading, each antenna 25 
stream is accumulated by an accumulator 240 and input into 
a bank of Fast Hadamard Transformers (FHTS) 250. Each 
FHT 250 produces N (e.g., N-64 for IS-95) complex values 
indicating the level of correlation between the values on that 
antenna stream and each of the N Walsh symbols. The 64 30 
M -element FHT output response vectors, y i7 are collected 
into an Mx64 matrix Y and processed by PSC 110 as in 
equation (1). Each response vector y i is cross-correlated with 
the array manifold data stored in local buffer memory 200. 
The number of complex multiply-and-accumulate (MAC) 35 
operations required per IS -95 Walsh symbol is 

Af*AT*K«16«64*256-262144, (5) 

assuming K-256 possible angles and M-16 antenna ele- 
ments. Since each IS -95 Walsh symbol has duration 208.33 40 
/isec (for a rate of 4800 Hz), PSC 110 must execute 
262144*4800=1.26*10° MAC operations per second. 

As described in above-referenced U.S. patent application 
Ser, No. 08/929,638 entitled "PRACTICAL SPACE-TIME 
RADIO METHOD FOR CDMA COMMUNICATION 45 
CAPACITY ENHANCEMENT", if a sufficient number (six 
or more) of antenna elements is used, the array manifold can 
be represented with a small number of bits. This reduced 
number of bits allows a simpler design of the PSC 110 to 
meet the 1.26* 10 9 MAC operations per second requirement 50 
and allows storage of the manifold coefficients in a smaller 
size memory 200. The result of the PSC correlation opera- 
tion is the KxN p matrix described above. 

The p matrix is further processed to produce a magnitude 
and angle of arrival (AO A) estimate each symbol time. From 55 
an implementation point of view, the most efficient process- 
ing on the P matrix is to simply select the overall maximum 
value. The corresponding row and column index provide the 
Walsh symbol (not needed when just doing AOA estimation) 
and AOA estimate, respectively, which are sent (along with 60 
the magnitude) to DSP unit 120 for processing. Instead of 
actually saving the P matrix, the maximum value can be 
saved from the stream of cross-correlation operations. 

In an alternative method of P matrix processing, the 
quality of the AOA estimate from PSC 110 for use in DSP 65 
unit 120 can be improved by averaging some matrix ele- 
ments before selecting the maximum value, During each 
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symbol time, a new K-element vector p is formed, whose 
elements are the maximum value from each row of matrix P. 
The vector p is then averaged with a number of p vectors 
from previous symbols, where the number of vectors to 
average is determined by the Doppler spread and is selected 
to correspond to a time period during which the spatial 
signature is relatively stable. During each symbol time, the 
maximum value of the new averaged vector is selected and 
sent, along with its position (AOA index), to DSP unit 120. 
Averaging before selecting a maximum value improves the 
signal-to-noise ratio on the AOA estimate. 

DSP unit 120 then creates and processes the AOA histo- 
gram and generates a matrix V, as described above. The 
matrix V is transferred to the EMP 130, which performs the 
operations described above and feeds the resulting beam- 
former coefficients to beamformer 140. Using the beam- 
former coefficient vector w, the beamformer 140 spatially 
filters the corresponding M-element vector y ( from the M 
FHT blocks 250. The N resulting values from the beam- 
former 140 give an indication of the likelihood of each of the 
N possible symbols. This information is then sent to the 
standard decoder blocks which are described in "CDMA, 
Principles of Spread Spectrum Communication" by A. J. 
Viterbi, referenced above. 

FIG. 3 shows a single finger of the coherent MAD suitable 
for 3G proposals (i.e., CDMA2000 and W-CDMA). A 
significant difference in the 3G proposals relative to the 
current IS -95 standard is the existence of a pilot signal 
embedded in the transmitted signal in the reverse link. 
Details of the reverse link structure for each proposal are 
given in "The cdma2000 ITU-R RTT Candidate Submis- 
sion" and "Japan's Proposal for Candidate Radio Transmis- 
sion Technology on IMT-2000: W-CDMA", respectively, 
incorporated by reference above. 

Similar to IS-95, each of M output signals from an 
antenna array are de-spread and accumulated in 
de-spreading blocks 300. However, unlike IS-95, each 
antenna stream is de-spread into two signals because the 3G 
systems incorporate an additional pilot signal along with the 
data signal. Even within the two 3G systems, de-spreading 
blocks 300 have slightly different structures due to the slight 
differences between the W-CDMA and CDMA2000 propos- 
als. 

FIG. 4 shows one of the M de-spreading blocks 300 for 
W-CDMA. The combined baseband signal is first 
de-scrambled by sending the signal to a complex multiplier 
410 for multiplication with an appropriate scrambling code. 
The de-scrambled signal is then sent to a set of second 
complex multipliers 420 for multiplication with a pilot 
spreading code to de-spread the pilot signal and with a data 
spreading code to de -spread the data signal, respectively, 
which separates the signal into respective pilot and data 
signals. The pilot and data signals are then transmitted to a 
set of accumulators 430, which accumulates the respective 
pilot and data signals over the number of pseudo-noise (PN) 
chips corresponding to one data symbol. Note that the pilot 
signal is rotated by 90° with respect to the data signal since 
the pilot is transmitted on the quadrature channel. 

FIG. 5 shows one of the M de-spreading blocks 300 for 
CDMA2000. The combined baseband signal is first input to 
a complex multiplier 510 for de-spreading with an appro- 
priate long code. The resulting signal is then transmitted to 
a second complex multiplier 520 for multiplication with an 
appropriate short code. The signal from second multiplier 
520 is sent to a third complex multiplier 530, which multi- 
plies the signal with the appropriate PN code to de -spread 
the signal and separate out the data signal. Note that the pilot 
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signal is spread with the all ones sequence. The pilot and 
data streams are sent to a pair of accumulators 540 for 
accumulation over the number of PN chips corresponding to 
one data symbol, resulting in received signal values corre- 
sponding to the pilot and data symbols. Note that in 5 
CDMA2000, the data signal is rotated by 90° with respect to 
the pilot signal since the data is transmitted on the quadra- 
ture channel. 

Referring back to FIG. 3, the M-element de-spread and 
accumulated received vectors y p and y rf , corresponding to 10 
the pilot and data signals, respectively, are transmitted to 
data estimator blocks 310 for the next step in the channel 
estimation process. The M data estimator blocks 310 use the 
y p and y d vectors to produce the two columns of the Mx2 
matrix Y For W-CDMA, the Y matrix is given by 15 

y-\yd+j'y P -yd+jy P l (6) 

where For CDMA2000, the Y matrix is given by 

Y=\y P +}y^ P -}yA- CO 20 

After forming the Y matrix, data estimator blocks 310 
send the Y matrix to PSC 110 for spatial correlation with the 
array manifold matrix C, as described above with respect to 
equation (1). 

The number of complex multiply- and-accumulate (MAC) 25 
operations required in PSC 110 per data symbol (only two 
possible data symbols in this case) is 

M*2*K=\ 6*2*256=8192, (8) 

assuming, as previously, K=256 possible angles and M=16 30 
antenna elements. For either W-CDMA or CDMA2000, the 
minimum data symbol duration is approximately 0.5 micro- 
seconds. Thus, the PSC must execute up to 8192*2* 10 6 = 
1.64*10 10 MAC operations per second. As described in the 
implementation for IS -95, if a sufficient number of antenna 35 
elements are used (e.g., M=6 or more), a low number of bits 
can be used to represent the array manifold, which allows a 
smaller and simpler PSC 110 to meet the 1.64*10 10 MAC 
operations per second requirement and allows storage of the 
manifold coefficients in a smaller size memory 200. The 40 
result of the correlation operation in PSC 110 is the P matrix 
having two columns and K rows. This Kx2 p matrix can now 
be processed as described above with respect to IS-95 (i.e., 
by selecting the overall maximum value or by averaging 
some matrix elements before selecting the maximum value) 45 
to produce results for DSP unit 120, 

Because the P matrix for 3G systems has only two 
columns, the complete P matrix can be stored more reason- 
ably than the larger P matrices for IS-95. Thus, the P matrix 
for 3G systems can also be processed by averaging all, 50 
instead of just some, of the P matrix elements before 
selecting the maximum value. During each symbol time, the 
P matrix is averaged with a number of P matrices from 
previous symbols, where the number of matrices to average 
is determined by the Doppler spread and is selected to 55 
correspond to a time period during which the spatial signa- 
ture is relatively stable. During each symbol time, the 
maximum value of the new averaged matrix is selected and 
sent, along with its position (AOA index), to DSP unit 120. 

As with IS-95 systems, DSP unit 120 then creates and 60 
processes the AOA histogram and generates a matrix ^, as 
described above. The matrix V is transferred to the EMP 
130, which performs the operations described above and 
feeds the resulting beainformer coefficients to beamfonner 
140. Using the beamformer coefficient vector w, the beam- 65 
former 140 spatially filters the M-element data vector y d 
from the de-spreading blocks 300. 
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Both the W-CDMA and CDMA2000 proposals recom- 
mend coherent demodulation. The beamformer coefficient 
vector w from EMP 130 is used directly to demodulate the 
de-spread data vector y d , i.e., the operation of the beam- 
former 140 is equivalent to coherent demodulation. The 
complex valued scalar quantity resulting from the demodu- 
lation operation, w^y^ is then sent to standard decoding 
blocks corresponding to the encoding and modulation 
blocks described in the W-CDMA and CDMA2000 propos- 
als. 

Forward link (FL) (i.e., base station to mobile station) 
channel estimation and demodulation are performed in a 
similar way as the reverse link (RL) using the extended array 
manifold method described above. The following descrip- 
tion illustrates a method and implementation for the 
CDMA2000 proposal. Performing channel estimation in the 
mobile station (MS) provides several benefits relative to 
performing channel estimation only in the base station (BS) 
on the reverse link and using this information for beam- 
forming on the forward link. Some benefits include provid- 
ing beamforming information when FL and RL usage is 
highly asymmetrical (i.e., when there are only infrequent RL 
transmissions), and providing more accurate beamforming 
information when angle spread is large, which can make it 
more difficult to determine the FL mean AOA from RL 
channel estimation. 

In the present invention, signals transmitted from each of 
the M antenna elements in the BS include a different 
auxiliary pilot spread with an extended Walsh code. This 
method is compatible with the CDMA2000 proposal. Use of 
dedicated pilots per MS is avoided since this is less efficient 
and leads to lower overall performance. The present inven- 
tion also allows the use of orthogonal transmit diversity 
(OTD), as described in the proposed CDMA2000 standard. 
Coded data bits in the BS are split into J data streams that 
are transmitted through corresponding J beams. Coefficients 
for the beams are determined using information from the MS 
as described below. 

FIG. 6 shows a single finger (demodulation channel) of a 
manifold assisted demodulator (MAD) 600 in the MS for the 
forward link of a CDMA2000 system. MAD 600 is used for 
both channel estimation and demodulation. Note that in 
contrast with the MAD finger for the reverse link, the EMP 
is used to filter the input data to the beamformers and not to 
filter the coefficients themselves. The received signal from 
the BS, which has been converted to baseband, is input to a 
de-spreading block 610 of the MS for separation of the pilot 
and data signals. 

FIG. 7 shows one embodiment of de-spreading block 610. 
The received combined signal is first transmitted to a 
complex multiplier 710 for de-spreading by complex mul- 
tiplication with the "PN spreading" or "short" code. The 
resulting signal is sent to a second complex multiplier 720 
for multiplication with the user "long" code, which 
de-scrambles the data signal. The output signal of second 
multiplier 720 is transmitted to a set of J complex multipliers 
730 for de -spreading the data signal, and the output signal of 
multiplier 710 is transmitted to a set of M complex multi- 
pliers 740 for de-spreading the pilot signal. 

The data signals are de-spread into J data streams by 
complex multiplication with the appropriate Walsh code. 
Thus, where OTD is being used, J different data streams are 
obtained. The effective number of spreading codes per user 
is the same as without OTD as a result of the splitting of the 
data streams and the extended Walsh code construction 
described in the proposed CDMA2000 standard. Each of the 
J de-spread data streams are then transmitted to a corre- 
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spo riding one of J accumulators 750 for accumulation over 
the number of PN chips corresponding to one data symbol. 

The orthogonal pilot signals are de-spread into signals 
corresponding to each of the M antenna elements by com- 
plex multiplication with corresponding extended Walsh 5 
codes. Each of the M de -spread pilot signals are transmitted 
to a corresponding one of M accumulators 760 for accumu- 
lation over a number of PN chips corresponding to a time 
period during which the spatial signature is relatively stable. 
In general, this time period is determined by the Doppler 10 
spread. The resulting M-element vector of pilot signal values 
is denoted as y fp . 

First, channel estimation techniques requiring only 
"slow" (i.e., <1 kbps) overhead messages between the BS 
and MS and vice versa are described. Identical beamformer 15 
coefficients are used in the BS and the MS, with switching 
of coefficients synchronized to signal frame boundaries. 

As part of the access procedure, the MS receives, from the 
BS, the array manifold table (with each coefficient repre- 
sented with a small number of bits as described previously) 20 
for that BS. The array manifold data is stored in a local 
buffer memory 620 and is used by a primary spatial corr- 
elator (PSC) unit 630 to perform the PSC operation 
expressed in equation (5), which is similar to equation (1), 
with Y and P matrices replaced with y fp and p^ vectors. 25 

pr&y* (9) 

The pf vector is processed as described above with the 
IS-95 system. The magnitude and angle of arrival (AO A) 
estimates each symbol time produced from PSC unit 630 are 30 
sent to a DSP unit 640 in the MS. DSP unit 640 generates the 
AO A histogram and determines the mean AOA value and 
angle spread estimate, similar to the processing described 
above. 

Once the AOA histogram has been generated and the 35 
mean AOA value and angle spread estimate have been 
determined, this information is used by DSP unit 640 to 
generate the subspace matrix V for EMP 650, as described 
above. This information is also used by DSP unit 640 to 
determine the beamforming coefficients, x, (i=l to J), for the 40 
J transmit beams. These beamforming coefficients x ( - will 
typically be selected as columns from the array manifold 
within the region of angle spread around the mean AOA 
value. These coefficients are sent to the BS to be used for 
transmit beamforming and are used by the MS in the 45 
beamformer blocks shown in FIG. 6. Note that in FIG. 6, the 
output signals from the EMP to the beamformer blocks are 
not the beamformer coefficients as described for the reverse 
link case. Instead, the EMP output signals are now the input 
signals to the beamformer blocks. The EMP output signals 50 
simply represent a filtered version of y^, which is now a 
better estimate of the forward link spatial signature. The 
filtered version of y^ is obtained from the output signals of 
the EMP, which is given by equation (4) (i.e., equation (4) 
where w represents the filtered version of y^,). 55 

Alternatively, the DSP unit can be divided between the 
MS and the BS. The MS generates the AOA histogram and 
associated mean AOA value and angle spread estimate and 
sends these values to the BS for further processing. The BS 
uses this information to generate the subspace for the EMP 60 
as described above. This information is also used to deter- 
mine the beamforming coefficients, x,- (i=l to J), for the J 
transmit beams. These beams will typically be selected as 
columns from the array manifold within the region of angle 
spread around the mean AOA value. These coefficients are 65 
used by the BS for transmit beamforming and are sent to the 
MS to be used in the beamformer blocks shown in FIG. 6, 
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In the case where "fast feedback" (e.g., 1 to 10 kbps rate) 
from the MS to the BS is feasible, actual spatial signature 
estimates (i.e., filtered versions of y fp ) can be used for 
beamforming coefficients in the BS instead of the coeffi- 
cients from the array manifold as described in the previous 
section. This is also possible when the forward link spatial 
signature is changing at a slower rate than would be typical 
for mobile applications (e.g., wireless local loop 
applications) such that the feedback rate becomes "slow" 
(<1 kbps). 

In this case, the filtered version of y fp at the output 
terminal of the EMP is sent back to the BS to be used as the 
coefficients for the transmit beamformer. In FIG, 6, for J-l, 
there is one beamformer block, and the coefficients of this 
one block are the same as the coefficients sent to the BS. The 
coefficients used in the BS and MS are synchronized so that 
they are changed at the same time. 

The next step in coherent demodulation of the data signal 
from each of the J beams is applying the beamformer 
coefficient vector x t - for beam i (i=l to J) to the spatial 
signature estimate of y fp from the EMP. This is just the inner 
product, c—x/y^, between the two vectors. In this case, 
however, the transpose of the coefficient vector is used, not 
the hermitian transpose as previously described. The result- 
ing complex scalar output signal, c^, is multiplied (after first 
conjugating) by the de-spread data streams, d,-, from 
de-spreading unit 610 to produce data estimates, c^d,, to be 
further processed in the standard CDMA2000 forward link 
decoder. 

Although the invention has been described with reference 
to particular embodiments, the description is only an 
example of the invention's application and should not be 
taken as a limitation. Various adaptations and combinations 
of features of the embodiments disclosed are within the 
scope of the invention as defined by the following claims. 

We claim: 

1. A method of signal processing for a wireless commu- 
nication system, comprising; 

receiving a set of de -spread vectors y„ wherein i=l to N 
and N is the number of possible received symbols from 
a wireless transmission; 

spatially correlating each of the de -spread vectors y ( with 
an array manifold matrix C to produce a correlation 
matrix P, wherein each vector y,- has M components, M 
being equal to the number of antenna output elements; 

determining angle of arrival (AOA) information of each 
vector y,- from matrix P; 

determining angular distribution information from a col- 
lection of AOA information; 

creating a subspace matrix from the matrix C and the 
AOA and angular distribution information; 

projecting each vector y, into the subspace matrix to 
produce a set of N vectors z,; 

finding a vector z^^ from the one of vectors z t - having the 
maximum energy; 

averaging vector z max to get a vector z avg ; and 

multiplying vector z avg with the subspace matrix to pro- 
duce a beamforming coefficient vector w. 

2. The method of claim 1, wherein the determining 
comprises finding the maximum element of the matrix P. 

3. The method of claim 1, wherein the determining 
comprises finding the maximum element of the matrix 
resulting from averaging a set of P matrices. 

4. The method of claim 1, further comprising, after the 
multiplying, determining the inner product of each vector y ( 
with the vector w. 
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5. The method of claim 1, wherein the array manifold determining the inner product, x^y fp7 of the vector y fp 
matrix C is an MxK matrix, K being equal to the number of with each beamforming coefficient vector x,- 

angles in the array manifold. 23. The method of claim 22, further comprising, after 

6. The method of claim 5, wherein the elements of the i th determining the inner product, multiplying the conjugate of 
column of the correlation matrix P represent a level of $ eacn mner product result with an associated one of J 
correlation between vector y, and the array manifold for each de-spread data signal values. 

of the K angles, 24. The method of claim 22, wherein pilot signals sent 

7. The method of claim 1, wherein creating the subspace from each antenna element in the base station are each 
matrix comprises a singular value decomposition of a por- spread with a different extended Walsh code. 

tion of the array manifold matrix C. 1Q 25. A beamforming system for wireless communication, 

8. The method of claim 1, wherein the subspace matrix is comprising: 

an approximation of the actual subspace. a primary spatial correlator (PSC), wherein the PSC 

9. The method of claim 8, wherein the approximation is receives de-spread data vectors y,- (i=l to N) and signals 
obtained by selecting columns of the actual subspace con- representing an array manifold matrix C and transmits 
taining dominant singular values. signals representing spatial information of the vectors 

10. The method of claim 1, wherein the averaging is over 15 y ; ; each vector y, having M components, M being equal 
a period approximately an order of magnitude shorter than to the number of antenna output elements and N being 
the Doppler spread of the received signals. equal to the number of possible symbols from a wire- 

11 . The method of claim 1, wherein the system is an IS-95 less transmission; 

based system. a digital signal processor (DSP) coupled to the PSC, 

12. The method of claim 11, wherein N is equal to 64. 20 wherein the DSP determines additional spatial infor- 

13. The method of claim 11, wherein the determining mation and determines a subspace containing most of 
comprises selecting the row and column of the matrix P the energy from the vectors y ; ; and 

containing the overall maximum value of the matrix p. an extended manifold processor (EMP) coupled to the 

14. The method of claim 11, wherein the determining D sp and the vectors y ; , wherein the EMP projects 
comprises: 25 signals representing each vector y ( - into the subspace to 

selecting the maximum-valued element from each row of generate projection coefficient vectors z i and outputs 

the matrix P to form a K-element p vector; signals representing a beamforming coefficient vector 

averaging the p vector; and , . ^ , • i . i • c 

. j- , 26. The system of claim 25, wherein the spatial informa- 

selecting the maximum value of the averaged p vector and ^ ^ ^ ang]< ? Qf anival (AQA) information . 

corresponding position. 27 ^ system of daim 2Sf wherein the additional spatial 

15. Hie method of claim 1, wherein the system is a information is the angular distribution information resulting 
W-CDMA or CDMA2000 proposed system. &om a collection of AOA information. 

16. The method of claim 15, wherein the y, vectors are a 2 8. The system of claim 25, wherein the PSC performs a 
combination of de-spread pilot vectors y p and de-spread data ^ spatial correlation between each of the vectors y { and the 
vectors y d . matrix C to produce signals representing a correlation matrix 

17. The method of claim 15, wherein N is equal to 2. p. 

18. The method of claim 16, wherein the combination 29. The system of claim 25, wherein the subspace matrix 
comprises a real part and an imaginary part. results from a singular value decomposition of a portion of 

19. The method of claim 16, wherein the determining the array manifold matrix C. 

comprises finding the maximum element of the matrix P. 40 30. The system of claim 25, wherein the subspace matrix 

20. The method of claim 16, wherein the determining is an approximation of the actual subspace. 

comprises finding the maximum element of the matrix 31- The s y stem of ciaim 30 > wherein the approximation 

resulting from averaging a set of P matrices. consists of columns of the actual subspace which contain 

21. The method of claim 16, wherein the determining do ™ n ™ { sin S ular values. m 

comprises averaging over all the elements of the matrix P 45 i 32 ". The s y stem ° f claim 25 > whercin the EMP comprises 

and finding the resulting maximum element and correspond- selectl0n circunrv f ™ electing a vector z_ having the 

• • maximum energy of the projection coefficient vectors z r 

m 22. Amethod of signal processing for a wireless comnra- »• ™ e s y stem of claim 31 < wherein the EMP further 



mcation system, comprising: 



comprises averaging circuitry for averaging vectors to 



i • , , •! , . ! 50 generate an average projection coefficient vector z flV „ that 

receiving, in a mobile station, a de-spread pilot signal & , . r , J * £ J , . , ^ 

, ' . . it _ ' . \ A r t ° u when multiplied by the subspace matrix produces an esn- 

vector y fr , wherein the vector ^ has M components, M 0 f the spatial signature 

being equal to the number of antenna output elements ^ ^ * 33 wherein ^ RMp 

in the base station; . w v e u- 1 • *u * <L 

, . . . . , comprises a multipher for multiplying the vector z with 

spatially correlating the vector y fp with an array manifold 55 ^ sub matrix t0 duce the vector w 

matrix C to produce a correlation vector P> 35 ^ system of daim ^ further comprising a beam . 

determining angle of arrival (AOA) information for the former coupled t0 the EMP and me signal vec tors y„ 

vector y/p from vector p/ wherein the beamformer determines the inner product, w^y,., 

determining angular distribution information from a col- 0 f eacn vector y ( with the vector w. 

lection of AOA information; 60 3$ The system of claim 25, wherein the system is an 

creating a subspace matrix from the matrix C and the IS-95 system, 

AOA and angular distribution information; 37. The system of claim 36, wherein N is equal to 64. 

projecting the vector y fp into the subspace matrix to 38. The system of claim 28, wherein the system is an 

produce a filtered version of the vector y fp ; IS-95 system, 

calculating beamforming coefficients x ( -, i-1 to J, from the 65 39. The system of claim 38, wherein the PSC selects the 

matrix C and the AOA and angular distribution infer- row and column of the matrix P containing the overall 

mation; and maximum value of the matrix P. 



11/01/2003, EAST Version: 1.4.1 



US 6,501,747 Bl 



17 



18 



40. The system of claim 38, wherein the PSC comprises: 
first selection circuitry for selecting the maximum-valued 

element from each row of the matrix P to form a 

K-element p vector; 
averaging circuitry for averaging the p vector; and 
second selection circuitry for selecting the maximum 

value of the averaged p vector and corresponding 

position. 

41. TTie system of claim 28, wherein the system is a 
W-CDMA or CDMA2000 proposed system. 

42. The system of claim 41, further comprising M data 
estimators, each data estimator receiving a de -spread pilot 
vector y p and de-spread data vector y d to form combinations 
of y p and y d to use as vectors y,-. 

43. The system of claim 42, wherein: 

the PSC further comprises a receiver for receiving a set of 
vectors y, from the data estimators and a transmitter for 
transmitting spatial information of the vectors y^; and 

the DSP is adapted to determine a subspace containing 
most of the energy from the vectors y,.. 

44. The system of claim 41, wherein N is equal to 2. 

45. The system of claim 41, further comprising a beam- 
former coupled to the EMP and the de -spread data vectors 
y d , wherein the beamformer determines the inner product, 
w^y^, of each vector y d with the vector w. 

46. A demodulation system for wireless communication 
comprising a plurality of manifold assisted demodulator 
(MAD) fingers, a first one of the plurality of MAD fingers 
comprising: 

a first block for de-spreading pilot and data signals and 
generating a de-spread pilot vector y p and a de-spread 
data vector y d \ 

a first data estimator block that forms complex valued 
combination vectors, y ( , of the de-spread pilot and data 
signal vectors; 

a first primary spatial correlator (PSC), wherein the PSC 
receives the de-spread data vectors y ( (i=l to N) and an 
array manifold matrix C and transmits spatial informa- 
tion of the vectors y„ each vector y, having M 
components, M being equal to the number of antenna 
output elements and N being equal to the number of 
possible symbols from a wireless transmission; 

a first extended manifold processor (EMP), wherein the 
first EMP receives the vectors y, and projects each 



vector y ( - into a subspace to generate projection coef- 
ficient vectors z,. and outputs a beamforming coefficient 
vector w; and 

s a first beamformer coupled to the first EMP and the data 
signal vectors y d , wherein the first beamformer deter- 
mines the inner product, w^y^, of each vector y d with 
the vector w. 

47. The system of claim 46, wherein a second one of the 
10 plurality of MAD fingers consists of: 

a second block for de -spreading the pilot and data signals 
and generating de-spread pilot vector y p and de -spread 
data vector y d ; 

15 a second data estimator block that forms the complex 
valued combination vectors, y iy of the de-spread pilot 
and data signal vectors; and 
a second primary spatial correlator (PSC), wherein the 
second PSC receives the de-spread data vectors y,- (i=l 
to N) and the array manifold matrix C and transmits 
spatial information of the vectors y t , each vector y i 
having M components, M being equal to the number of 
antenna output elements and N being equal to the 
number of possible symbols from a wireless transmis- 
sion. 

48. The system of claim 47, wherein a third one of the 
plurality of MAD fingers consists of: 

a third block for de-spreading the pilot and data signals 
and generating de-spread pilot vector y p and de-spread 
data vector y d ; 
a third data estimator block that forms the complex valued 
combination vectors, y,-, of the de -spread pilot and data 
signal vectors; and 
a second extended manifold processor (EMP), wherein 
the second EMP receives the vectors y, and projects 
each vector y ( into a subspace to generate projection 
coefficient vectors z t and outputs a beamforming coef- 
ficient vector w; and 
a second beamformer coupled to the second EMP and the 
data signal vectors y di wherein the first beamformer 
determines the inner product, w"y d , of each vector y d 
with the vector w. 
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